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magnetoresistance: spin-orbit coupling in magnetic tunnel junctions 
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The effects of the spin-orbit interaction on the tunneling magnetoresistance of ferromag- 
net/semiconductor/normal metal tunnel junctions are investigated. Analytical expressions for the 
tunneling anisotropic magnetoresistance (TAMR) are derived within an approximation in which the 
dependence of the magnetoresistance on the magnetization orientation in the ferromagnet originates 
from the interference between Bychkov-Rashba and Dresselhaus spin-orbit couplings that appear at 
junction interfaces and in the tunneling region. We also investigate the transport properties of ferro- 
magnet/semiconductor/ferromagnet tunnel junctions and show that in such structures the spin-orbit 
interaction leads not only to the TAMR effect but also to the anisotropy of the conventional tun- 
neling magnetoresistance (TMR). The resulting anisotropic tunneling magnetoresistance (ATMR) 
depends on the absolute magnetization directions in the ferromagnets. Within the proposed model, 
depending on the magnetization directions in the ferromagnets, the interplay of Bychkov-Rashba 
and Dresselhaus spin-orbit couplings produces differences between the rates of transmitted and re- 
flected spins at the ferromagnet/seminconductor interfaces, which results in an anisotropic local 
density of states at the Fermi surface and in the TAMR and ATMR effects. Model calculations for 
Fe/GaAs/Fe tunnel junctions are presented. Furthermore, based on rather general symmetry con- 
siderations, we deduce the form of the magnetoresistance dependence on the absolute orientations 
of the magnetizations in the ferromagnets. 

PACS numbers: 73.43.Jn, 72.25.Dc, 73.43.Qt 
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I. INTRODUCTION 



The tunneling magnetoresistance (TMR) effect is 
observed in fcrromagnct/insulator/ferromagnet hetcro- 
junctions, in which the magnetoresistance exhibits a 
strong dependence on the relative magnetization direc- 
tions in the two ferromagnetic layers and on their spin 
polarizationsF^^^ Because of this peculiarly strong 
asymmetric behavior of the magnetoresistance, TMR de- 
vices find multiple uses ranging from magnetic sensors to 
magnetic random access memory applications^^ 

Beyond the conventional TMR effect, it has been 
observed that the magnetoresistance in magnetic tun- 
nel junctions (MTJs) may also depend on the abso- 
lute orientation of the magnetizations in the ferromag- 
netic leadsi&i&S This phenomenon is called the tunneling 
anisotropic magnetoresistance (TAMR) effect ££ A the- 
oretical investigation of the tunneling magnetoresistance 
in GaMnAs/GaAlAs/GaMnAs tunnel junctions was re- 
ported in Ref. J9J which predicted that, as a result of 
the strong spin-orbit interaction the tunneling magne- 
toresistance depends on the angle between the current 
flow direction and the orientation of the electrode magne- 
tization. A difference between the tunneling magnetore- 
sistanccs in the in-plane (i.e., magnetization in the plane 
of the magnetic layers) and out-of-plane configurations of 
up to 6% was predicted for large values of the electrode 
spin polarization^ Here we refer to this phenomenon as 
the out-of-plane TAMR. Recent first principles calcula- 
tions in Fe/MgO/Fe magnetic tunnel junctions (MTJs) 
predict an out-of-plane TAMR ratio of about 44 %J^ On 
the other hand, we refer to an in-plane TAMR effect as 



the change in the magnetoresistance when the in-plane 
magnetization of the ferromagnetic layer(s) is rotated in 
the plane perpendicular to the direction of the current 
flow. 

It is remarkable that the TAMR is present even in 
MTJs in which only one of the electrodes is magnetic and 
the conventional TMR is absent. In contrast to the con- 
ventional TMR-bascd devices, which require two mag- 
netic layers for their operation, TAMR-based devices can 
operate with a single magnetic lead, opening new possi- 
bilities and functionalities for the operation of spintronic 
devices. The TAMR may also affect the spin-injection 
from a ferromagnet into a non-magnetic semiconductor. 
Therefore, in order to correctly interpret the results of 
spin injection experiments in a spin-valve configuration, 
it is essential to understand the nature, properties, and 
origin of the TAMR effect. 

The first experimental observation of TAMR 
was in (Ga,Mn)As/A10x/Au hetcrojunctions, in 
which an in-plane TAMR ratio of about 2.7% was 
founds Experimental investigations of the in-plane 
TAMR in (Ga,Mn)As/GaAs/(Ga,Mn)As and in 
(Ga,Mn)As/ZnSe/(Ga,Mn)As tunnel junctions, in 
which both electrodes are ferromagnetic have also been 
reported^ In the case of (Ga,Mn)As/ZnSe/(Ga,Mn)As, 
the in-plane TAMR ratio was found to decrease with 
increasing temperature, from about 10% at 2 K to 
8.5% at 20 K<£ This temperature dependence of 
the in-plane TAMR is more dramatic in the case of 
(Ga,Mn)As/GaAs/(Ga,Mn)As, for which a TAMR 
ratio of the order of a few hundred percent at 4 K was 
amplified to 150 000% at 1.7 Ki 7 - This huge amplification 
of the in-plane TAMR was suggested to originate from 



the opening of the Efros-Shklovskii gap 1 - 1 - at the Fermi 
energy when crossing the metal-insulator transition. 7 
Measurements of the TAMR in p+-(Ga,Mn)As/n+- 
GaAs Esaki diode devices have also been reported^ 
In addition to the investigations involving vertical 
tunneling devices the TAMR has also been studied 
in break junctions ) 13 ! 14 nanoconstriction o 12 ' 15 and 
nanocontactsJ^ 

Beyond the area of currently low Curie temper- 
ature ferromagnetic semiconductors, the TAMR 
has recently been experimentally investigated in 
Fc(001)/vacuum/bcc-Cu(001) tunnel junctions^ 
Fe/GaAs/Au MTJs^ Co/AlO^/Au MTJsA CoPt 
structures^ and in multilayer-(Co/Pt)/A10 x /Pt 
structures^ 

In what follows we focus our discussion on the case 
of the in-plane TAMR (for brevity we will refer to 
it as the TAMR effect). We investigate the TAMR 
in ferromagnct/scmiconductor/normal metal (F/S/NM) 
and in fcrromagnet/semiconductor/ferromagnet (F/S/F) 
MTJs. We propose a model in which the two-fold sym- 
metric magnctorcsistancc dependence on the orienta- 
tion of the in-plane magnetization in the ferromagnetic 
layer(s) originates from the interference of Dresselhaus 
and Bychkov-Rashba-like spin-orbit couplings. Such in- 
terference effects have already been investigated in lat- 
eral transport in 2D electron systcm a 22 i 23 ' 24 , in spin re- 
laxation in quantum wells 2 ^ and quantum dots 2 ^, or in 
2D plasmonsi^ 7 . The symmetry, which is imprinted in the 
tunneling probability becomes apparent when a magnetic 
moment is present. Our main results are: i) finding an- 
alytical expressions for evaluating the TAMR in both 
F/S/NM and F/S/F MTJs, ii) prediction and evaluation 
of the ATMR in F/S/F heterojunctions, and iii) deriva- 
tion of a simple phcnomenological relation describing the 
dependence of the tunneling magnetoresistance on the 
absolute orientation of the magnetization(s) of the ferro- 
magnet(s). 

The paper is organized as follows. In Scc.[Tl]we present 
the theoretical model describing the tunneling through a 
MTJ. In a first approximation we consider the case of 
an infinitesimally thin barrier (Sec. lII A)l . while the finite 
spatial extension of the potential barrier is incorporated 
in a more sophisticated approach discussed in Sec. Ill Bl 
Detailed solutions and tunneling properties within these 
approximations are given in Appendices |X] and [B] rcspec- 
tivel y. In S ecHTTlwc discuss the TA MR in both F/S/NM 
(Sec. HE! and F/S/F (Sec. EH) MTJs. The ATMR in 
F/S/F tunnel junctions is investigated in Sec. HVi where 
specific calculations for model Fc/GaAs/Fe MTJs are 
presented. In Sec. [V] we develop a phenological model for 
determining the dependence of the TAMR and ATMR on 
the absolute oricntation(s) of the magnctization(s) in the 
ferromagnetic lead(s). Finally, conclusions are given in 
Sec. El 



II. THEORETICAL MODEL. 

Consider a F/S/F tunnel heterojunction. The semi- 
conductor is assumed to lack bulk inversion symmetry 
(zinc-blende semiconductors are typical examples). The 
bulk inversion asymmetry of the semiconductor together 
with the structure inversion asymmetry (for the case 
of asymmetric junctions) of the heterojunction give rise 
to the Dresselhaus5i22i22j2S and Bychkov-Rashba^^ML 
SOIs, respectively. The interference of these two spin- 
orbit couplings leads to a net anisotropic SOI with a C<i v 
symmetry which is imprinted onto the tunneling mag- 
netoresistance as the electrons pass through the semi- 
conductor barrier. This was discussed in some details 
in Refs. [alia for the case of F/S/NM tunnel junctions. 
Here we generalize the model proposed in Refs. [5lfl8l to 
the case of F/S/F tunnel junctions. For such structures 
our model predicts the coexistence of both the TAMR 
and ATMR phenomena. 

We consider a F/S/F tunnel junction grown in the 
z = [001] direction, where the semiconductor forms 
a barrier of width d between the left and right ferro- 
magnetic electrodes. At first we discuss a simplified 
model for very thin barriers. In that case the bar- 
rier can be approximated by a Dirac delta function and 
the SOI reduced to the plane of the barrier. In what 
follows we will refer to this model as the Dirac delta 
model (DDM). A second model in which Slonczewski's 
proposal 3 ^ for ferromagnet/insulator/ferromagnet tun- 
nel junctions is generalized to the case of ferromag- 
net/semiconductor/fcrromagnct junctions by including 
the Bychkov-Rashba and Dresselhaus SOIs will be re- 
ferred to as the Slonczewski spin-orbit model (SSOM). 



A. Dirac delta model (DDM) 

We consider here the case of a very thin tunneling bar- 
rier. Assuming that the in-plane wave vector k|| is con- 
served throughout the heterostructure, one can decouple 
the motion along the growth direction (z) from the other 
spatial degrees of freedom. The effective model Hamilto- 
nian describing the tunneling across the heterojunction 
reads 



H = H + H z + H l 



so- 



Here 
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2m,Q dz 2 



V dS(z), 
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(2) 



with too the bare electron mass, and Vq and d the high 
and width, respectively, of the actual potential barrier 
[here modelled with a Dirac delta function S(z)] along 
the growth direction (z = [001]) of the heterostructure. 

The spin splitting due to the exchange field in the left 
(z < 0) and right (z > 0) ferromagnetic regions is given 



by 



H z = - 



©(- 



-n, 



Q(z)A r 
a n r 



(3) 



Here Aj and A r represent the exchange energy in the 
left and right ferromagnets, respectively, and ®(z) is the 
Heaviside step function. The components of the vector 
cr arc the Pauli matrices, and \\j = (cos#j, sin#j,0) with 
j = l,r is a unit vector defining the in-plane magneti- 
zation direction in the left (j = I) and right (j = r) 
ferromagnets with respect to the [100] crystallographic 
direction. The Zeeman splitting in the semiconductor 
can be neglected. 

We note that in recent experiments with Fc/GaAs/Au 
tunnel junctions^ the reference axis was taken as the 
[110] direction. Therefore, it is convenient to express 
the magnetization direction relative to the [110] axis by 
introducing the angle shifting <fij = 9j — ir/4 (j = l,r). 
One can then write rr, = [cos(<pj + 7r/4), sin((j>j + 7r/4), 0] 
with <pj giving the magnetization direction in the left 
(j = I) and right (j = r) ferromagnets with respect to 
the [110] crystallographic direction. 

Within the DDM, the spin-orbit interaction through- 
out the semiconductor barrier (including the interfaces) 
can be written as 

H so = (w • <r)5{z), (4) 

with the effective spin-orbit coupling field 

w = {-ak y + ^/k x ,dk x - 7fc y ,0). (5) 

Here a and 7 represent effective values of the Bychkov- 
Rashba and linearized Dressclhaus parameters, respec- 
tively, and k x and k y refer to the x and y components 
of the wave vector k. In terms of the usual Dresselhaus 
parameter 7, the linearized Dresselhaus parameter can 
be approximated^ as 7 « 7Q, where Q = 2moVod/H 2 
stands for the strength of the effective wave vector in the 
barrier. 

The scattering states in the left (z < 0) and right (z > 
0) ferromagnetic regions are given by 
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respectively. Here we have introduced the wave vector 
components 




and 




(8) 



(9) 



with fc|| = w fc^ -(- fc2 denoting the length of the wave 

vector component corresponding to the free motion in 
the x — y plane. The spinors 
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tt/4) U = hr), 



(10) 



correspond to a spin parallel (a =|) or antiparallcl 
(cr =1) to the magnetization direction rij = [cos(</>j + 
7r/4), sm(cj)j + 7r/4) , 0] in the left (j = I) and right (j = r) 
ferromagnets. 

The reflection and transmission coefficients can be 
found by imposing appropriate boundary conditions and 
solving the corresponding system of linear equations (for 
details see Appendix \Q . The transmissivity of an in- 
coming spin-er particle can then be evaluated from the 
relation 

TviE,^) =Rc[K a (\t a , a \ 2 + K „ a \t a ^\ 2 )}. (11) 

Explicit analytical expressions for the transmission coef- 
ficients (icr.CT and t a ^ a ) are given in Appendix IA1 



B. Slonczewski spin-orbit model (SSOM) 

We give a generalization of the Slonczewski model^ 
for ferromagnct/insulator/fcrromagnet tunnel junctions 
to the case in which the insulator barrier is replaced by 
a zinc-blende semiconductor. Unlike in the DDM, now 
the spatial extension of the potential barrier is taken into 
account. The model Hamiltonian is 



H = H a + H z + H 



BR 



Hr 



where 



#o 



m(z) 



V e(z)Q(d-z). 



(12) 



(13) 



The electron effective mass m(z) is assumed to be m = 
m c in the central (semiconductor) region and m = ma 
in the ferromagnets. The exchange splitting in the ferro- 
magnets is now given by 



e(-z)A, e(z - d)A r 

H z = n; • a- n r • a. 

The Dressclhaus SOI can be written as^ ' 29 ' 30 ' 32 ^ 



(14) 



d 
H D = (k x a x - kyGy)— ( 7(2) 



d_ 

' dz 



(15) 



where x and y correspond to the [100] and [010] direc- 
tions, respectively. The Dressclhaus parameter "f(z) has a 
finite value 7 in the semiconductor region, where the bulk 
inversion asymmetry is present, and vanishes elsewhere. 
Note that because of the step-like spatial dependence of 
7(2), the Dresselhaus SOI [Eq. ([15)) ] implicitly includes 
both the interface and bulk contributionsi 5 ' 29 i 33 



The Bychkov-Rashba SOI is given by^s 



H B r = [aiS(z - zi) - a r 5{z - z r )\ (k x a v - k y a x ), 



(16) 



and arises due to the ferromagnct/semiconductor inter- 
face inversion asymmetry^ Here a; (a r ) denotes the SOI 
strength at the left (right) interface z; = (z r = d). 
For the small voltages considered here (up to a hundred 
mV) , the Bychkov-Rashba SOI inside the semiconductor 
can be neglected. 

The z-components of the scattering states in the left 
and right ferromagnets have the same form as in Eqs. (J6j> 
and {7J, respectively 

In the central (semiconductor) region (0 < z < d) we 
have 



¥^ = Y J (A^^ Z + B 



- 9 * z )xf\ 
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The angle £ is defined through the relation tan(£) = 
—hy/kx- We have also used the notation 



<1\ 



where 



*) 



Qa 






/it(^ 


■y 


l2m c (V - 


E) 


+ q 


h? 





(19) 



(20) 



is the length of the z component of the wave vector in 
the barrier in the absence of SOI. 

The expansion coefficients in Eqs. |6]), ([7]), and (fT7l can 
be found by applying appropriate matching conditions at 
each interface and by solving the corresponding system of 
linear equations (for details, see Appendix [B]) . Once the 
wave function is determined, the particle transmissivity 
can be calculated from Eq. (jlip . Approximate analyti- 
cal expressions for the transmission coefficients t a , a and 
ta.-a are given in Appendix IBl 



III. TUNNELING ANISOTROPIC 
MAGNETORESISTANCE (TAMR) 

The magnctorcsistance of a tunnel junction can be ob- 
tained by evaluating the current through the device or 
the conductance. The current flowing along the hetero- 
j unction is given by 



I 



(2TT) 3 h 



J2 I 'd£d 2 k||T CT (£,k|| )[/,(£)- f r (E)], 

(21) 
where fi(E) and f r (E) are Fermi-Dirac distributions 
with chemical potentials fii and /i r in the left and right 
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FIG. 1: (color online), (a) Schematics of a Fe/GaAs/Au MTJ. 
The magnetization direction in the ferromagnet is specified 
by the vector n. (b) Schematics of the potential profile of the 
heterojunction along the [001] direction, (c) Definition of the 
TAMR in F/S/NM junctions. 



(metallic or ferromagnetic) leads, respectively. For the 
case of zero temperature and small voltages, the Fermi- 
Dirac distributions can be expanded in powers of the 
voltage Vbias ■ To first order in Vuas one obtains /; (E) — 
f T {E) « 5(E — Ep)Vbias with S(x) the Dirac delta func- 
tion and Ep the Fermi energy. One then obtains the 
following approximate expression for the conductance 



G 



o-=T,J. 



Go- , G CT 



(2tt) 3 ^ 



y"d 2 k||T CT (£; F ,k||). (22) 



We note that although similar, the expression above 
differs from the linear response conductance. In our 
case, the transmissivity T a {Ep,')s.n) depends on the 
Bychkov-Rashba parameter a. Recent first-principles 
calculations^ have shown that the spin-orbit coupling 
field is different for different bands, therefore the effec- 
tive value of a is energy dependent. By applying an ex- 
ternal voltage the energy window relevant for tunneling 
can be changed, resulting in voltage-dependent values of 
a. Consequently, the conductance in Eq. (|22)l depends, 
paramctrically, on the applied voltage. 



A. TAMR in ferromagnet/semiconductor/normal 
metal tunnel junctions 

The tunneling properties of F/S/NM junctions can be 
obtained as a limit case of the models proposed in Sec. 
|TT]for F/S/F tunnel junctions by taking 0; = <f> r = <f> and 
A r as the Zeeman splitting in the normal metal region. 
In the present case I, c, and r refer to the ferromagnetic 
(left), semiconductor (central), and normal metal (right) 
regions, respectively. 

The TAMR in F/S/NM tunnel junctions refers to the 
changes of the tunneling magnctorcsistance (R) when 



varying the magnetization direction n; of the magnetic 
layer with respect to a fixed axis. Here we assume the 
[110] crystallographic direction as the reference axis. The 
TAMR is then given by£ 



TAMR [11O] (0) 



R(4>) - R(0) _ G(0) - G(4>) 



R(p) 



G{<£) 



(23) 



Since in a F/S/NM tunnel junction only one electrode is 
magnetic, the conventional TMR effect is absent. 

An alternative to the magnetoresistance, which refers 
to the charge transport, is the spin polarization efficiency 
of the transmission characterized by the tunneling spin 
polarization^ 



P = 



h-h 



(24) 



where I a is the charge current corresponding to the spin- 
cr channel and / is the total current. The changes in 
the tunneling spin polarization when the magnetization 
of the ferromagnet is rotated in-plane can then be char- 
acterized by the tunneling anisotropic spin polarization 
(TASP), which is defined as£ 



TASP [110] O) 



P(0) ~ P(<t>) 



(25) 



Taking into account that the Zeeman splitting in 
the normal metal is small we can approximate K a s=s 
K-a- Then for the DDM the conductance is given by 
Eq. (TA"T9|) . It follows from Eqs. fl23|) and Eq. (TAl9|) that 
the TAMR is given by 



TAMR [llo] (0) 



G s 



>(0) - G £ 



\<f>) 



G iso + G amso(^) 



(26) 



For junctions in which the Bychkov-Rashba and Dressel- 
haus SOIs can be considered as small perturbations, the 
anisotropy is small and G amso ((/>) <C /G 1S0 . In addition, 
the spin-orbit contribution G' = Gi + G 1 ^ [see Eqs. (|A8[) - 
(|A10[) ]^1 to the isotropic part of the conductance is also 



much smaller than the contribution G^ 



g; o) 



G 



(0) 



corresponding to the system in the absence of the SOI, 
i.e., G' <C G(°>. Therefore one can drop the contributions 
G' and G amso from the denominator in Eq. $M. The 
TAMR can then be approximated as 



TAMR [110] (</>) 



G 8 



'(0) - G £ 



\4>) 



G(°) 



(27) 



The substitution of Eq. (|A17[) into Eq. (gT]) leads to 
TAMR[ 110] (^)w 

,2 



S2 T fcjj) + (ff21^ 



G(°) 



A Q A 7 [cos(20)-l]. (28) 



Here we have introduced the dimcnsionlcss SOC param- 
eters A a = 2mQa/h 2 and A 7 = 2mo7//i 2 . The functions 
g2\ and gi\ are given by Eq. (|A18|) . 



[-110] 




FIG. 2: (color online), (a) Schematics of the anisotropy of 
the spin-orbit magnetic field w(ky). Thin arrows represent 
a vector plot of the SOC field w. The solid line is a polar 
plot of the SOC field strength |w(kii)| for a fixed value of 
fey = |k|||. When the magnetization of the ferromagnet points 
along the [-110] direction [see thick blue (dark) arrow], the 
direction of the strongest SOC field is parallel to the incident 
majority spins which easily tunnel through the barrier. On 
the contrary, for a magnetization direction [110] [see thick 
green (light) arrow], the strongest SOC field is perpendicular 
to the incident spins and the tunneling becomes less favorable. 
The net result is a spin valve effect whose efficiency depends 
upon the absolute orientation of the magnetization and gives 
rise to the TAMR. 



The expression above gives the angular dependence of 
the TAMR and is consistent with the angular depen- 
dence experimentally observed in Fa/GaAs/Au tunnel 
junctions^ It also suggests that bias-induced changes 
of the sign of the Bychkov-Rashba parameter a result 
in the inversion (change of sign) of the TAMR (such 
an inversion has been experimentally observed^) Fur- 
thermore, one can see from Eq. (|2"8]) that the amplitude 
of the TAMR is governed by the product A a A 7 ex 0:7 
and the averages (g2ak?,)cr {a =1\4)- When 0:7 = 0, 
the two-fold TAMR is suppressed (the suppression of the 
TAMR was also observed in Reffla). i.e., as long as other 
anisotropic effects such as uniaxial strain are not present, 
the Bychkov-Rashba (or Dresselhaus) SOI alone cannot 
explain the experimentally observed Civ symmetry of the 
TAMR. The TAMR vanishes also if the spin polarization 
of both electrodes becomes sufficiently small. In such a 
case £;f,t ~ ^f,x and gi a vanishes [sec Eq. (|A18|) ]. result- 
ing in the suppression of the TAMR. On the contrary, 
Eq. (f2"5| predicts an increase of the TAMR amplitude 
for F/S/NM tunnel junctions whose constituents exhibit 
large values of «7 as well as a large spin polarization in 
the magnetic electrode. 

A simple, intuitive explanation of the origin of the uni- 
axial anisotropy of the TAMR can be obtained by investi- 
gating the dependence of the effective spin-orbit coupling 
field w(k||) [see Eq. (JSJ)], i.e., the effective magnetic field 
that the spins feel when traversing the semiconducting 



barrier. A schematics of the anisotropy of the spin-orbit 
field w(k||) is shown in Fig. (J5J), where the thin arrows 
represent a vector plot of w(k||), while the solid line is a 
polar plot of the field amplitude |w(k|| )| for a fixed value 
of fen = |k|||. The spin-orbit field is oriented in the [110] 
([-110]) direction at the points of low (high) spin-orbit 
field, where the field amplitude |w| reaches a minimum 
(maximum). When the magnetization in the ferromag- 
net points along the [-110] direction, the direction of the 
highest spin-orbit field is parallel to the incident, major- 
ity spins which are then easily transmitted through the 
barrier. On the other hand, for a magnetization direction 
[110], the highest spin-orbit field is perpendicular to the 
incident spins (to both the majority and minority spins) 
and the transmission becomes less favorable than in the 
case the magnetization is in the [-110] direction. This 
spin-orbit induced difference in the tunneling transmis- 
sivities depending on the magnetization direction results 
in the uniaxial anisotropy of the TAMRj2& Furthermore, 
the magnetization direction dependence of the transmis- 
sion and reflection of the incident spins should be re- 
flected in the local density of states at the interfaces of 
the barrier. Within the DDM the left (F/S) and right 
(S/NM) interfaces are merged into a single plane and one 
can not distinguish between them. A more detailed view 
of the role of the interfaces requires the use of the SSOM. 
It turns out (this will be shown latter in this section) that 
the F/S interface plays a major role in the TAMR phe- 
nomenon while the S/NM interface appears irrelevant. 
This is intuitively expected since the exchange splitting 
in the ferromagnet is much larger than the Zeeman split- 
ting in the normal metal. Consequently, the spin-valve 
effect at the F/S interface is much stronger than in the 
S/NM interface. 

The local density of states reflect also the uniaxial 
anisotropy of the TAMR with respect to the magneti- 
zation orientation in the ferromagnet. In fact, one can 
introduce the anisotropic, local density of states (AL- 
DOS) through the definition 

Ainnq ( a\ LDOS(z,O)-LDOS(z,0) 

ALDOS[no](2;,<?>) = tt^c^.. j,\ > ( 29 ) 



LDOS(z,« 



where 



LDOS(z,0)= £ |^L |* ff (*, 



<j>,k„ F )\ , (30) 



"=t,i 



is the total local density of states at position z and evalu- 
ated at the Fermi surface determined by the Fermi wave 
vectors 



k a F — 



'2m 



E, 



4 1 



(31) 



Since we are interested only in propagating states, we 
may restrict the possible values of k\\ to the interval 



[0, kf nax ], with k° nax given by Eq. (|A16|) . Since the spin 
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FIG. 3: (color online). Values of the ratio -Rrno] /-Rrno] as a 
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splitting in the normal metal region is negligibly small, 
K a ss K_ a . It follows from Eqs. Q and (fTTj) that 



T CT (£,fc||)oc|a>S r) | 2 , 



(32) 



and, therefore, the conductance [see Eq. (|22p] is related 
to the LDOS at z = d as G a {<f) ex LDOS(d, (j>). One then 
obtains that 



TAMR [11O] (0) « ALDOS [110] (z = d,<f>). 



(33) 



For a numerical illustration we consider an epitaxial 
Fe/GaAs/Au heterojunction similar to that used in the 
experimental observations reported in Ref. [l8|. We use 
the value m c = 0.067 tuq for the electron effective mass 
in the central (GaAs) region. The barrier width and 
hight (measured from the Fermi energy) are, respectively, 
d = 80 A and V c = 0.75 eV, corresponding to the ex- 
perimental samples in Ref. Qj3. For the Fe layer a Stoner 
model with the majority and minority spin channels hav- 
ing Fermi momenta fcp-f = 1.05 x 10 8 cm -1 and kpi = 
0.44 x 10 8 cm _1 j21 respectively, is assumed. The Fermi 
momentum in Au is taken as up = 1.2 x 10 8 cm" 1 ^ 
We consider the case of relatively weak magnetic fields 
(specifically, B = 0.5 T). At high magnetic fields, say, 
several Tesla, our model is invalid as it does not in- 
clude cyclotron effects which become relevant when the 
cyclotron radius approaches the barrier width. 

The Dresselhaus spin-orbit parameter in GaAs is 7 ~ 
24 eV A ££&£. On the other hand, the values of the in- 
terface Bychkov-Rashba parameters ai , a r [sec Eq. (flrJj) ] 
arc not know for metal-semiconductor interfaces. Due 
to the complexity of the problem, a theoretical estima- 
tion of such parameters requires first-principles calcula- 
tions including the band structure details of the involved 
materials^ which is beyond the scope of the present pa- 
per. Here we assume cti and a r as phenomenological pa- 
rameters which must be understood as the values of the 
interface Bychkov-Rashba parameters at the ferromag- 
net/semiconductor and semiconductor/normal metal in- 
terfaces, respectively, averaged over all the relevant bands 
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FIG. 4: (color online). Angular dependence of the TAMR in 
a Fe/GaAs/Au MTJ. Solid and dashed curves represent the 
results obtained within the SSOM and the DDM, respectively, 
(a) for voltages -90 and -50 mV. (b) for voltages 50, 90, and 
135 mV. 



contributing to the transport across the corresponding in- 
terfaces. In order to investigate how does the degree of 
anisotropy depend on these two parameters we performed 
calculations of the ratio R\iio] / R[uo] (which is a measure 
of the degree of anisotropyi^) as a function of cti and a r 
by using the spin-orbit Slonczewski model described in 
Sec. Ill Bl The results are shown in Fig.[3j where one can 
appreciate that the size of this ratio (and, consequently, 
of the TAMR) is dominated by ai . This is because the 
Zeeman splitting in Au is very small compared to the 
exchange splitting in Fe and, consequently, the spin flips 
mainly when crossing the ferromagnet/semiconductor in- 
terface. Then, since the values of the TAMR are not very 
sensitive to the changes of a r we can set this parameter, 
without loss of generality, to zero. This leaves ai as a 
single fitting parameter when comparing to experiment. 
The values of the phenomenological parameter at were 
determined in Refs. [5|[l8| by fitting the theory to the ex- 
perimental value of the ratio -fir-no] /-Rfiio] and a very 
satisfactory agreement between theory and experiment 
was achieved. This fitting at a single angle was enough 
for the theoretical model to reproduce the complete an- 
gular dependence of the TAMR, demonstrating the ro- 
bustness of the model. 

Assuming that the interface Bychkov-Rashba param- 
eter ai is voltage dependent and performing the fitting 
procedure for different values of the bias voltage the bias 
dependence of ai can be extracted^^ Here we use the 
same values of ai reported in Refs. [afl8l for computing 
the angular dependence of the TAMR at different val- 
ues of the bias voltage. The results are shown in Fig. 5J 
where the dashed and solid lines correspond to calcula- 
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FIG. 5: (color online). Angular dependence of the TASP in 
a Fe/GaAs/Au MTJ for different values of the bias voltage. 



tions within the DDM and the SSOM, respectively An 
overall agreement between the two models can be appre- 
ciated. The TAMR exhibits an oscillatory behavior as a 
function of the magnetization direction [see also Eq. (J2HJ)] 
and can be inverted by changing the bias voltage [com- 
pare Figs. |4{a) and (b)]. This bias induced inversion of 
the TAMR was experimentally observed in Fe/GaAs/Au 
tunnel junctions^ and was explained to occur as a con- 
sequence of a bias induced change in the sign of the ef- 
fective Bychkov-Rashba parameter. Preliminary ab ini- 
tio calculations^ for Fe/GaAs structures suggest that 
the Bychkov-Rashba parameters associated with differ- 
ent bands may have different values and even change the 
sign. Thus, the effective value of the interface Bychkov- 
Rashba parameter cti will depend on which bands are the 
ones that mainly contribute to the transport across the 
Fe/GaAs interface (at low temperature those are the ones 
which have the appropriate symmetry and lie inside the 
voltage window around the Fermi energy). For different 
values of the bias voltage different set of bands will be 
relevant to transport. On the other hand, to different set 
of bands correspond different effective values of the in- 
terface Bychkov-Rashba parameter. Consequently, a; be- 
comes strongly dependent on the bias voltage. The above 
analysis leads to the conclusion that the inversion of the 
TAMR originates from the bias induced sign change of 
the effective Bychkov-Rashba SOI at the Fe/GaAs inter- 
face, as proposed in Refs. [afl8l. 

The sign change of the effective Bychkov-Rashba pa- 
rameter has also influence on the tunneling spin polar- 
ization, resulting in the bias induced inversion (change 
of sign) of the TASP as shown in Fig. [SJ Bias in- 
duced changes of the sign of the tunneling spin polar- 
ization in Fc/GaAs/Cu MTJs has also been reported^ 
The anisotropy of the tunneling spin polarization, which 
also exhibits a two-fold symmetry, indicates that the 
amount of transmitted and reflected spin at the inter- 
faces depends on the magnetization direction in the Fe 
layer, resulting in an anisotropic local density of states 
at the Fermi surface. This is consistent with previous 
work a 6 i 7 i 21 in which the anisotropy of the density of states 
with respect to the magnetization direction was related 
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FIG. 6: (color online), (a) Angular dependence of the TAMR 
in a Fe/GaAs/Au MTJ for different values of the barrier width 
d and a bias voltage V bias — —90 mV. (b) Amplitude of the 
TAMR (at = 90°) as a function of the barrier width. The 
TAMR curves were obtained by using the SSOM. 



to the origin of the TAMR. In fact, our model calcula- 
tions reveal that the TAMR, ALDOS [see Eq. (|33|)]. and 
TASP all exhibit a C2 V symmetry with the same kind of 
angular dependence (compare Figs. [4] and [5]). 

A system parameter that can influence the size of the 
TAMR is the width of the barrier. The angular depen- 
dence of the TAMR calculated within the SSOM for the 
case of Vbias = —90 meV is displayed in Fig. GDJa) for dif- 
ferent values of the barrier width d. As one would expect, 
the changes in the barrier width do not affect the two-fold 
symmetry of the TAMR but only its amplitude, whose 
absolute value is predicted to increase when increasing 
the width of the barrier [see Fig. [5Jb)]. 

We remark that our model neglects the contribution of 
the spin-orbit-induced symmetries of the involved bulk 
structures. Say, Fc exhibits a four-fold anisotropy, which 
should be reflected in the tunneling density of states. The 
fact that this is not seen in the experiment suggests that 
this effect is smaller than the two-fold symmetry consid- 
ered in our model. 



B. TAMR in 

ferromagnet /semiconductor /ferromagnet tunnel 
junctions 



Our discussion in Sec. IIII Al suggests that in the case 
of a F/S/F tunnel junction, the magnetoresistance will 
depend on the absolute direction of the magnetization 
in each of the ferromagncts. In the left and right elec- 
trodes, the magnetization directions with respect to the 
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[110] crystallographic direction are given by the angles 
4> — 4>i and <j) r , respectively. For convenience we intro- 
duce the angle 9 = 4> r — <pi, describing the magnetiza- 
tion direction in the right ferromagnet relative to that in 
the left ferromagnet (see Fig. [7]). Different values of the 
tunneling magnetoresistance are expected to occur when 
in-plane rotating the magnetizations of both ferromag- 
nets at the same time, while keeping the relative angle 9 
fixed. Thus, the expression for the TAMR in F/S/NM 
junctions [see Eq. ([25)) ] can now be generalized as 



TAMR [110] ( 



R{B,4>)-R{0,Q) G(O,O)-G(0,t 



R(6,0) 



G(M) 



(34) 

Following the same procedure as in Sec. IIII Al one can, 
in principle, obtain analytical expressions for the TAMR 
in a F/S/F junction. It turns out however that in the 
general case defined by Eq. (f3"4")) the resulting relations 
are quite lengthy and not much simpler than the more 
accurate expressions obtained within the SSOM (see Ap- 
pendix [B]). Therefore, we omit here the expressions re- 
sulting from the DDM and show only the results obtained 
within the SSOM. 

The dependence of the TAMR on the angles 9 and </> is 
shown in Fig. g] for the case of a Fe/GaAs/Fc MTJ. The 
width of the barrier is d = 80 A. Two different cases, 
corresponding to a > (a) and a < (b), arc consid- 
ered. In both cases, the absolute value of the TAMR 
reaches its maximal amplitude when the magnetization 
of the left electrode is parallel to the [-110] direction (i.e., 
<j) = 90°, 270°) and the one of the right electrode is per- 
pendicular to it (i.e., 9 = 90°, 270°). This is because this 
configuration, in our parabolic model, corresponds to the 
case of the stronger structure inversion asymmetry and, 
consequently, to the largest absolute value of the effec- 
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FIG. 8: (color online). Calculated TAMR denned in Eq. d34j) 
in a Fc/GaAs/Fe MTJ as a function of the angles <j> and 9. 
(a) for ai = 42 eV A 2 and a r = 32 eV A 2 (i.e., a > 0). (b)for 
m = 42 eV A 2 and a r = 52 eV A 2 (i.e., a < 0). 



tive Bychkov-Rashba parameter a. We also note that at 
a fixed value of 9 the TAMR has a two-fold symmetry 
with respect to 4> and vice versa. It is clear, however, 
that the angles 9 and <f> play different roles in the sym- 
metry of the TAMR, which manifest in the lack of mirror 
symmetry with respect to the axis 9 — <$>. We have in- 
vestigated some traces of the TAMR displayed in Fig. [8] 
The results are shown in Fig. [5] where we present polar 
plots of the TAMR as a function of <j> for fixed values of 
9. The solid lines correspond to the calculations within 
the SSOM. The meaning of the dashed lines will be ex- 
plained in Sec. [V] The orientation of the symmetry axis 
of the two-fold symmetric TAMR is determined by the 
relative angle 9 rather than by the relative values of the 
interface Bychkov-Rashba parameters (note for the same 
9, the symmetry axis does not change its orientation with 
a r ). The amplitude of the TAMR is bigger for the case 
of a r = 52 eV A 2 than for a r = 32 eV A 2 and shows a 
strong dependence on 9. 



IV. ANISOTROPIC TUNNELING 
MAGNETORESISTANCE (ATMR) 

We now consider the SOI effects on the TMR for 
the case of a F/S/F tunnel junction. The conven- 
tional TMR effect in ferromagnct/insulator/fcrromagnct 
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FIG. 9: (color online). Calculated angular dependence of the 
TAMR in a Fe/GaAs/Fe MTJ for fixed values of the angle 
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(F/I/F) tunnel junctions relies on the dependence of 
the magnetoresistancc across the junction on the rela- 
tive magnetization directions in the different ferromag- 
netic layers and their spin polarizations^ and is usually 
defined as 



TMR 



Rap — Rp 
R~p ' 



(35) 



where Rp (Rap) is the magnctorcsistance measured 
when the magnetization of the left and right ferromag- 
netic layers are parallel (antiparallcl). 

In the case of F/S/F heterojunctions, the interference 
between Bychkov-Rashba and Dressclhaus SOIs leads to 
anisotropic effects. Consequently, the conventional TMR 
in F/S/F junctions will depend not only on the relative 
but also on the absolute magnetization directions in the 
ferromagnets, resulting in an anisotropic magnctorcsis- 
tance. For quantifying this phenomenon we use the fol- 
lowing anisotropic generalization of the TMR, 



ATMR [110] = 



Rap&) - Rp(<I>) Gp(4>) - G A p{ct>) 



Rp{4>) 



Gap{4>) 



(36) 

which now accounts for the magnetoresistancc depen- 
dence on the absolute magnetization orientations with 
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respect to the [110] crystallographic direction (see Fig. [7]). 
Furthermore, the efficiency r\ of the anisotropic effects on 
the tunneling magnetoresistance can be defined as 



ATMR [11O] (0)-ATMR [11O] (O) 
ATMR [110] (0) 



(37) 



A simplified approximate expression for the ATMR can 
be found within the DDM by following the same proce- 
dure as in Sec. IIII Al The result is 



ATMR [110] 



TMR 



hG 



(0) 

AP 



^4-^E ((9^-9t P )kf 



°-=U 



A Q A 7 cos(20) 



Z^ 

°=l,l 



i9£-&,W\ 



(38) 



where TMR 



{Gf 



,(0) 



,(0) 



G A p)/G Ap is the conventional 



92a, 



TMR in the absence of SOI, and the functions g^ 
gff, and g^f are given by the Eqs. l[A"26 ]) . (|X27 j) . (jA"33|) . 
and (|A34[) . respectively. Assuming that the spin-orbit 
effects on the conventional TMR are small one can ap- 
proximate the ATMR efficiency as 



A a A 7 [cos(20) — 1] 



hG 



(o) 

AP 



TMR 



E 

'=1,1 



(92a 



9l P )k\ 



(39) 

Note that the angular dependence of the efficiency r\ is 
similar to that of the TAMR given in Eq. (f28|) . 

It follows from Eq. (f3"B")) that when both Bychkov- 
Rashba and Dresselhaus SOIs are present (i.e., when 



A Q A 7 



a l 7^ 0), the magnetoresistance becomes 



anisotropic, resulting in the ATMR effect here predicted. 
One can see also that this effect exhibits a two-fold sym- 
metry with a cos(2</>) angular dependence. Unlike the 
TAMR, when aq = the ATMR becomes isotropic but 
docs not vanish. However, in such a limit the efficiency 
?/ of the ATMR do vanish [see Eq. (JHSJ)] , indicating the 
absence of anisotropy. Like for the TAMR, changes of 
the sign of a result in the inversion (change of sign) of 
the efficiency 77 of the ATMR. On the other hand, since 
the TMR contribution in Eq. (j3"5)) is usually the dom- 
inant one, the bias-induced changes of the sign of the 
Bychkov-Rashba parameter a will not, in general, cause 
the inversion (change of sign) of the ATMR. However, 
the fact that the TMR contribution may change sign in 
dependence of the applied voltag o 39 i 40 ' 41 can result in the 
inversion of the ATMR. 

We performed calculations within the SSOM of the an- 
gular dependence of the ATMR for a Fe/GaAs/Fe MTJ 
with a fixed value of the left-interface Bychkov-Rashba 
parameter (a; = 42 eV A 2 ) corresponding to a bias volt- 
age V bia ^ = —90 mV . The results are shown in Figs. [TO] 
(a) and (b) for the parameters a r = 32 eV A 2 and 
a r = 52 eV A 2 . The angular dependence of the ATMR 
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FIG. 10: (color online). Angular dependence of the ATMR 
defined in Eq. ((36} for a Fe/GaAs/Fe MTJ. (a) for ai = 42 
eV A 2 and a r = 32 eV A 2 (i.e., a > 0). (b)for a, = 42 eV A 2 
and a r = 52 eV A 2 (i.e., a < 0). 



is consistent with Eq. (|3"8|) . The orientation of the sym- 
metry axis of the ATMR is determined by the sign of 
the effective value a. On the other hand, the amplitude 
of the ATMR shows a behavior opposite to the one cor- 
responding to the TAMR (see Fig. \§§, in the sense that 
now the amplitude of the anisotropy of the TMR is big- 
ger for the case a r = 32 eV A 2 than for a r = 52 eV A 2 
and, unlike for the case of the TAMR (see Fig. [5]), the 
orientation of the symmetry axis of the ATMR can be 
flipped by varying the value of a r . 



V. PHENOMENOLOGICAL MODEL OF TAMR 
AND ATMR 

In order to explain the origin of the angular depen- 
dence of the TAMR [sec Eq. (25)] in F/S/NM junctions, a 
phenomcnological model based on rat her g eneral symme- 
try considerations wasproposed in RefQia and elaborated 
in more details in Ref|5J. Here we extend this phenomeno- 
logical model to the case of F/S/F tunnel junctions. 

For a given k|| , there are three preferential directions 
in the system: i) the magnetization direction n; in the 
left ferromagnet, ii) the magnetization direction n r in the 
right ferromagnet, and iii) the direction of the effective 
spin-orbit field w. A scalar quantity such as the trans- 
missivity for the a spin channel can be expanded in a 
series of the all possible scalars one can form with the 
vectors n/, n r , and w. Thus, to the second order in |w|, 
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we have 



T a = Ti 0) + T {2) . 



(40) 



Note that because of symmetry reasons,— the linear in |w| 
terms vanishes after averaging over k\\ and, therefore, we 
have omitted them in Eq. (|40|) . The zero order terms 
have the general form 



T(°)=4 j + 4 J(n z .„ r ), 



(41) 



,(0) 



.,«>) 



where a\J and a^J are isotropic expansion coefficients. 

Taking into account that n; = [cos(c/) + 7r/4),sin(0 + 
tt/4),0] and n r = [cos(0 + 4> + 7r/4),sin(0 + 4> + 7r/4),0], 
the general angular dependence of T^ can be extracted 
from Eq. (glj. The result is 



r<°> = 4°j + 4°J ooB(fl) 



(42) 



This equation describes the dependence of the transmis- 
sivity on the relative angle 9 between the magnetizations 
in the left and right ferromagnetic electrodes in the ab- 
sence of SOI and is consistent with previous results 4£ 

The second order contribution T^ can be cast in the 
following general form 



T (2 



[«2 

4 2) K 



(2; 



(n, ■ n r 



42 

6(7 



(n; • w) 



,(2) 



■w)|"+4; J ( n/ -w)(n r -w), (43) 

(2) 

with the expansion coefficients a ia (i = 1,2, ...,5) being 
angular independent. To second order in the SOI, the 
conductance is then given by 



G a — — 
h 



(<Ti°)) k|| + {Tf\ 



(44) 



where (...)k,, denotes average over kn and evaluation at 
the Fermi energy Ep. 

Taking into account Eqs. ©, ([42]). (|43|) . and (J44J) and 
performing the corresponding angular averaging in It- 
space one obtains the following, general form of the total 
G\ + Gi, 



conductance G 
G{6,4>) = 

where 



Gi + G 2 cos(6>) H 
G 4 cos(20 + 26») 



G 3 cos(20) 
-G 5 cos(20 + 6>), (45) 



and 



Gi = -r E (afJ^xWykl 



G 5 = "^ E \ at 5a W xWykf\ 



o=T,l 



(49) 



(50) 



In the equations above the average (...)<r have the same 
meaning as in Eq. (|A15|) . We note that for systems which 
are isotropic in the absence of the SOI these equations 
are quite general. They are valid up to second order in 
the strength of the spin-orbit coupling field regardless of 
the specific form of w. 

One can see that for the case of a F/S/NM junction or 
for the case of parallel and antiparallcl configurations in a 
F/S/F junction, the general angular dependence given in 
Eq. (|45p is consistent with the corresponding expressions 
obtained in Appendix [X] We want to stress, however, 
that the relations presented in Appendix [Al are the result 
of a specific approximation (the DDM), while Eqs. ([45]) 
- (1501) arc general, and valid for any model or approx- 
imation (as far as the system without SOI is isotropic, 
the effective spin-orbit field is weak, and n;, n r , and w 
lie in a plane perpendicular to the current flow). More 
general relations corresponding to arbitrary orientations 
of ii/ , n r , and w will be given elsewhere 42 

The general form of the TAMR in F/S/NM junctions 
follows from Eqs. ([5]), (T2"3"|) and (|4"5)) . Proceeding as in 
Sec. HE1 one finds 



TAMR [11O] (0)« 


e 2 1 


E (-(42 


h G(°) 


A Q A 7 [cos 


(20) - 1], 



,(2) 

'-1<T 



Ha ) K \\ 



(51) 



which is consistent with the corresponding expression 
found within the DDM [see Eq. ((28)) ]. Similarly, one ob- 
tains for the TMR and ATMR in F/S/F junctions the 
following general expressions, 



TMR 



1 



E 



h a (0) 



2n (0) k 2 



(52) 



Gi 



and 



E 



,(0) 



-Oi 



(2) 



(2) 



°-=T,i 



7 (2) ^ 2 
l 4a ) K \\ 



G2 -tE(4°. ) + ^(24 2 ct ) + 4 2 ^ 



(46) 



(47) 



G* = — 



e v--\ 



h <^ 



(o&WyAfi) , (48) 



ATMR [11O] (0) = 

e 2 (A 2 +A 2 ) 



TMR 



2h 



G 



(o) 



E(( 4 4 2 . ) + 24 2 . ) )fc 2 



h fifC) 



AP a=1,i 



,( 2 )l.2 



E ( 2 <>n 



respectively. Note that Eq. 
previous result in Eq. 



A Q A 7 cos(2</>), (53) 



is consistent with our 
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The TAMR in a F/S/F exhibits a more complicated 
angular dependence, which has the general form 



TAMR [m (B,<t>) = 

A(6)[l - cos(2</>)] + B(9) sin(20) 
Gi + G 2 cos + A{9) cos(20) - B{9) sin(20) 



where 



and 



A{9) = G 3 + G 4 cos(20) + G 5 cos 9, 



B(9) = G 4 sin(20) + G 5 siii0. 



,(54) 



(55) 



(56) 



At the phcnomenological level, the expansion constants 
Gi, ...G5 in Eqs. (j45|) and ([54]) can be determined from 
conductance measurements (or conductance theoretical 
evaluation) at selected values of 9 and </>. For example, 
the values of the phcnomenological constants can be de- 
termine as follows 



Gi = 



G(tt/2, 0) + G(tt/2, tt/2) 



G 2 = G(0,tt/4)-Gi 



G 5 = Gi-G(7r/2,7r/4), 



(57) 



(58) 



(59) 



In order to check the validity of the general angular 
dependence given in Eqs. (|4"5|) and ([54]) we have com- 
puted the expansion coefficients from Eqs. ([ST)) - ([6"Tj) by 
using the SSOM. We then evaluate the TAMR by us- 
ing the phenomenological relation given in Eq. (|54[) and 
compare the results (see dashed lines in Fig. O with the 
full angular dependence obtained from the SSOM (solid 
lines in Fig. \§§ . The agreement is satisfactory although 
some small discrepancies appear. We attribute these dis- 
crepancies to the fact that when determining the coeffi- 
cients Gi, ..., G5 from Eqs. ([STj) - ([6"Tj) we have ignored, 
as discussed above, that these coefficients are weakly 0- 
dependent. 

VI. CONCLUSIONS 



We have investigated the spin-orbit induced anisotropy 
of the tunneling magnetoresistance of F/S/NM and 
F/S/F MTJs. By performing model calculations we have 
shown that the two-fold symmetry of the TAMR effect 
may arise from the interplay of Bychkov-Rashba and 
Dresselhaus SOIs. The spin-orbit interference effects in 
F/S/F MTJs lead to the anisotropy of the conventional 
TMR. Thus, the ATMR depends on the absolute mag- 
netization directions in the fcrromagncts. The magne- 
toresistance dependence on the absolute orientation of 
the magnetization in the ferromagnets is deduced from 
general symmetry considerations. 



G 3 = G(tt/4, 0) - Gx 



V2 



and 



(G 2 + G 5 ) 



G 4 = G 1 + G 3 -G(7r/2,0). 



(60) 



(61) 



It is worth noting that although we have referred to 
the coefficients Gi, ..., G5 as constants, strictly speaking, 
these parameters may exhibit a dependence on the rela- 
tive angle 0. Such a dependence may originate from the 
fact that the averages containing the components of the 
spin-orbit field [see Eqs. (|4*6")) - ([5"0[) ] are, in principle, 0- 
dependent. However, as shown below, this effect is weak 
for the system here considered. 



J 



Acknowledgments 

We are grateful to M. Gmitra for useful hints and dis- 
cussions and to D. Weiss and J. Moser for discussions on 
TAMR related experiments. This work was supported 
by the Deutsche Forschungsgemeinschaft (DFG) through 
the Sonderforschungsbereich (SFB) 689. 



APPENDIX A 

The cigenfunctions of the Hamiltonian given by Eq. ([T]) 
obey the following matching conditions 



*«(0-) = v4 r) (0+); 



h 2 rf* 



(0 



2mo dz 



+ (Vod + w-<x)*«(0 



z=0" 



h 2 cti> 



(r) 



2rriQ dz 



(Al) 



z=0+ 



r 



which, with the scattering states in Eqs.Q and ([7|), lead coefficients are found to be given by 
to a system of 4 linear equations for determining the co- 



efficients 



cr.er , 1 a, — a s 



t CT ,cr; and t a -<j- The transmission 



4d 2 v^ (fe_ CT + K-a + %Q) (1 + e*<*-*i>) 



n 



+ 



and 



myr a (u ■ s g|g ) 



Ad\fk a ~d (k- a + K a + iQ) (1 - 



(A2) 



n 



8i (U ■ S n 



O 



where 



with 



n = n + (-)n_(+) - n+(+)n_(-), 



(A3) 



(A4) 



n±(A) = d(k ±a + K Xa + iQ)(l±Xe i ^-'l" 
+ 2z(U-S ±a , Aff ). 

The vectors So- jCr ' (a' = ±<j) are given by 



W<n&\ 



(A5) 



(A6) 



while Q = 2m Vod/h 2 and U = (2m rf/S 2 )w. 

The transmissivity for an incident particle with spin a 
is given by Eq. (fTTj) . For the case of a F/S/NM junction 
we can approximate k = n a ~ W— cr and the transmissiv- 
ity reduces to 



T a = Rc[K(\t a . a \ 2 + \t a .^\ 2 )] 



(A7) 



In order to obtain a simplified analytical expression for 
the TAMR in F/S/NM junctions, we consider the case 
in which the effective spin-orbit field is small. In such a 
case one can expand Eq. (|A7[) in powers of w = |w| and 
obtain the conductance from Eq. (|2"2"|) . The result, up to 
second order in w, reads 

G CT « GT + G^ niso . (A8) 

The isotropic part of the conductance is given by 

GT = G^ + G'^ (A9) 

where 



G (o) = l _ 
h 



4\k F ^\n F 



13 



with 



4 ±(T = fcF,±a + K F + iQ, (All) 

is the conductance in absence of the SOI and 



n' — " f\2 



(A£ + ^)<ffia*f>a, 



(A12) 



is the isotropic contribution induced by the SOI. Here we 
have denoted k Fa = k a (E F ,k\\), k f = K(E F ,kn) and 

8K F \k F , a \ ( 2Bc[A a A- a \ \ 

^-WlXZ^K l^l 2 - 1 )- 92 - (A13) 

We have also introduced dimensionless SOC parameters 
as 



2m _ 2m _ 

Ar, = . „ a, A~, = ^„ 7- 



ri 2 



h 2 



(A14) 



The function 92 is defined below in Eq. (|A18j) . The av- 
erage 



1 fKl 

1 / max 



2Wo 



k\\dk 



II' 



(A15) 



where 



fc max = min 



'2m / A, 



2?rio 



£/■ 



(A16) 
denotes the maximum values of k» for which we have inci- 
dent and transmitted propagating states. In the average 
defined in Eq. (|A15|) the corresponding angular integra- 
tion over ky [see Eq. (|2"2"|) ] has already been performed. 

The spin-orbit induced anisotropy of the conductance 
is determined by the relation 



G 



amso 



{4>) = --rKA 1 {g2akfj a cos{2(l)), 



(A17) 



(A10) where 



92a 



AK F \k F ^\ (\A_ a \ 2 \A a - A_ a \ 2 - 4Im[A_ ff ](|A_ ff | 2 Im[A CT ] - \A a \ 2 Im[A_ a })) 



I 



(A18) 



The total conductance can then be written as 

G(0) « G iso + G aniso (<^), (A19) 



with G iso = Gf ° + Gf° and G aniso (0) = Gf" so (0) 
G aniso (0). 
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We now consider the case of a F/S/F tunnel junction. 
For the general case of arbitrary orientations of the mag- 
netization in the left and right ferromagncts, the ana- 
lytical expressions for the conductance are very lengthy 
and we therefore omit them here. Simpler expressions 
are found, however, for the particular cases of parallel 
(P) and antiparallcl (AP) configurations. If we consider 
that the ferromagnetic electrodes are made of the same 
material, then k a = K a and the transmissivity reduces to 



and 






2h 



The anisotropic contribution reads 



(A24) 



T a = Re [k a \t a ^\ 2 + k- \t v ,- \ 2 ] 



(A20) 



Following the same procedure as above, the following ex- 
pression for the conductance in the parallel configuration 
is found 



Cr ib °W = X A Q A 7 ^J (g£kp a cos(2<t>). 

<T=T,i 



(A25) 



The isotropic part is given by 



(A21) 



The functions g^ and g^ a are given by 



^-.iso _ f~i 

\J p — Crp 



(0) 



G', 



where 



o-=T4 



pi 



4P 



4*2> 



( A22 ) P _ 8fep> (2k F , a + k F ,- a ) Q 2 - 32ftg > fa?,- a . p 

9la ~ ; Z2~/ L T 9la 



(A23) and 



4^,, + Q 2 Uk*+Q 2 



(A26) 



P _ 4fcf> (fcf.-g - fc F;Cr ) r 



A- 3 



12k Fi<T (k F ,a + k F ^ a )Q 2 



(A27) 



respectively. 

On the other hand, for the case of antiparallel config- 
uration one obtains the following expressions 



where 



Gai 



amso 



AP "+" ^AP 



(4>), 



,^iso _ ,o(o) , ,^,/ 

"AP — U AP "+" u APi 



(A28) 



(A29) 



and the functions g^ and g^f are given by 

Ap _ 8fcF,ff {2k F - a + k F . a ) Q 2 AP 
9io 



and 



(7 


\k F ,„ + k F ^ a f + Q 2 


3 y^a 


8k F ,a (k F ,cr ~ 2fcp_ (T ) (k.F,a + k F - a ) 



(fcp jCr + kF-a) + 



, (A33) 



fi(o) _ £_ y^ / 4fc F;g . (fcj^ - fc F;Cr ) 
" h a tt d \(k F .-, + k F .,) 2 + Q 2 



'AP 



(A30) 



92a 



Ak F ^ a (fcp CT - fe^.-o-) 



(*J 



-.) 2 + 



(A34) 



and 



G^p = -^(Xl + A 7 ) J2 (9i P k 2 )a- (A31) 



The anisotropic contribution for the antiparallel case 
reads 



GSp 1S0 (</>) = X A Q A 7 £ (gt p k\) a cos(2</>), (A32) 



respectively. 

In the derivation of Eqs. (|A24j) . (|A25|) . (|A31|) . and 
()A32[) we assumed that the spin orbit parameters a and 
7 are independent of the relative magnetization orienta- 
tions in the ferromagncts. Strictly speaking the values 
of a and 7 are ^-dependent and therefore to the parallel 
and antiparallel configurations correspond slightly differ- 
ent sets of effective spin orbit parameters. This effect 
turns out to be negligible for the system here investi- 
gated and we therefore omit it. 
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In both the parallel and antiparallcl cases, the averages 
have the same meaning as in Eq. (|A15[) but now with 



k° v defined as 



'2777-0 



n 2 \ EF + a Y 



(A35) 



Fe/GaAs/Fe structures, for example, one of the Fc/Gas 
interface may be Ga terminated while the other may be 
As terminated . 43 ' 44 This interface-induced asymmetry is 
enough for the average Bychkov-Rashba parameter to 
have a sizable value. 



It is worth noting that for a junction with structure 
inversion symmetry the average Bychkov-Rashba param- 
eter vanishes. One may think that since we are consid- 
ering a structure with leads of the same material the 
anisotropic term in the conductance corresponding to 
the parallel configuration (for the parallel configuration 
the structure becomes symmetric under spatial inversion) 
will vanishes. This is so for strictly symmetric under 
spatial inversion structures. In practice, however, the in- 
terfaces may not be identical, e.g., one of the interfaces 
may be epitaxial while the other not. Another possibil- 
ity is to consider different terminations of the semicon- 
ductor at the left and right interfaces. For the case of 



J 



APPENDIX B 

Here we discuss the details about the solutions of the 
SSOM and provide approximate expressions for the tun- 
neling coefficients. 

By requiring the probability flux conservation across 
the interfaces one obtains that the eigenfunctions of the 
Hamiltonian in Eq. (|T2|) should fulfil the following bound- 
ary conditions 



*£W =*£>(*«), 



(Bl) 



h 2 

2mj 



2m, 7i 



yrCx&x ky^y) 



'i.i 



(k x a y - k y a x )^'(zij), 



(M 



W 



dz 



2m i 



1- 



2m jlj 



\Kx@x ky&y) 



dVZ 



GO 



dz 



(B2) 



r 



where i,j = l,c,r and z; c = and z cr = d are the lo- 
cations of the left and right interfaces, respectively. The 
interface Bychkov-Rashba parameters are introduced as 
a ic = a i, a cr = -«n while the Dresselhaus parameter is 
7 C = 7 in the semiconductor and vanishes elsewhere (i.e., 

ll = Jr = 0). 

Applying the above boundary conditions to the scat- 
tering states given in Section IIIBI one obtains a system 
of 8 linear equations for determining the 8 unknown ex- 
pansion coefficients. The exact expressions for the trans- 



J 



mission coefficients t aM and t a _ a are quite cumbersome. 
However, simplified analytical expressions for t aM and 
£<7 _<7 arc found in the limit q$d 3> 1. In such a case one 
finds the following approximate relations for the tunnel- 
ing coefficients, 



tn 



D n 



t, 



D ' ">'" D ' 

where D = /_(-)/+(-) - /-(+)/+(+), with 



(B3) 



/±(A) = ±- 



2% 



( ae *(*r-ew4) + Ae 2iA 



2 \ m ±A 
I 



q TX - lK± A(l-Xae^^+^) 



(B4) 



and 

1 J_ / 2m c 7fc|| 

m± m c \ h 2 

Furthermore, we have 



(B5) 



D„ 



2mod 

m + 



q-f-(+)9- 



2mod 



and 



2m d 2m d 

q-f+{-)g- q+f+{+)g + . B7) 

m + m_ 



q+f-(—)g+, (B6) In these equations we introduced the notation 
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idJk„ 



9± 



(/o T hi - ^<z T ) (1 ± CTe *(^+«+-/4)) T fta (i T ae <( 



i !+ £+7r/4)> 



-g ± d 



fc2+(/o-fci 



mod 



<i- 



fo + hi-^ 



r 



'7+ 



(B8) 



where /o = i(/c CT + k- a )d/2, 



I'l 



and 



iad 



(k a - k_ a ) cos((j>i + £ + tt/4) + 



h-2 = -^r(K 



fc_ CT )sin(<^+£ + 7r/4). 



aik\\ 



Vo 



■sin(2£), 
(B9) 



-» " cos(2g). 

^0 



(BIO) 

It is worth noting that the approximate expressions 
for the tunneling coefficients here provided are valid up 
to first order in exp(— q±d). This approximation is ap- 
propriate for treating junctions with high and not too 
thin potential barriers. For the systems here considered 



the hight of the barrier (with respect to the Fermi level) 
is about Vb = (Vq — Ep) ss 0.75 cV and d varies from 
20 A to 200 A. In such cases the approximations here 
discussed turns out to be excellent. 

It is not difficult to show that in the limit ai = a r = 
7 = 0, the expressions for the tunneling coefficients here 
obtained reduce to the ones reported in Ref. for the 
case of vanishing SOI. 

We also remark that the expressions above were ob- 
tained for the general case of a F/S/F tunnel junction 
but the corresponding expressions for a F/S/NM junction 
can easily be obtained by taking the limits n a = K- a = k 
and d>i = d> r = 4>- 
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